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We investigate a U(l) lattice chiral gauge theory (LxGT) with domain wall fermions and gauge fixing. In the 
reduced model limit, our perturbative and numerical investigations at Yukawa coupling y — 1 show that there are 
no extra mirror chiral modes. The longitudinal gauge degrees of freedom have no effect on the free domain wall 
fermion spectrum consisting of opposite chiral modes at the domain wall and the anti-domain wall which have 
an exponentially damped overlap. Our numerical investigation at small Yukawa couplings {y <C 1) also leads to 
similar conclusions as above. 



1. Introduction 

Lattice regularization of chiral gauge theories 
has remained a long standing problem of nonper- 
turbative investigation of quantum field theory. 
Lack of chiral gauge invariance in LxGT propos- 
als is responsible for the longitudinal gauge de- 
grees of freedom (dof) coupling to fermionic dof 
and eventually spoiling the chiral nature of the 
theory. The obvious remedy to control the lon- 
gitudinal gauge dof is to gauge fix with a target 
theory in mind ||l|. The formal problem is that 
for compact gauge-fixing a BRST-invariant par- 
tition function as well as (unnormalized) expec- 
tation values of BRST invariant operators van- 
ish as a consequence of lattice Gribov copies [D . 
Shamir and Golterman |0] have proposed to keep 
the gauge-fixing part of the action BRST non- 
invariant and tune counterterms to recover BRST 
in the continuum. In their formalism, the contin- 
uum limit is to be taken from within the broken 
ferromagnetic (FM) phase approaching another 
broken phase which is called ferromagnetic direc- 
tional (FMD) phase, with the mass of the gauge 
field vanishing at the FM-FMD transition. This 
was tried out in a U(l) Smit-Swift model and so 
far the results show that in the pure gauge sector, 
QED is recovered in the continuum limit |^ and 
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in the reduced model limit free chiral fermions in 
the appropriate chiral representation are obtained 

i- 

When one gauge transforms a gauge non- 
invariant theory, one picks up the longitudinal 
gauge degrees of freedom (radially frozen scalars) 
explicitly in the action. The reduced model is 
then obtained by making the lattice gauge field 
unity for all links, i.e., by switching off the trans- 
verse gauge coupling. The reduced model is ob- 
viously a Yukawa model. The job of gauge fixing 
in the LxGT when translated into the reduced 
model is to find a continuous phase transition 
where these unwanted scalars are decoupled leav- 
ing only free fermions in the appropriate chiral 
representation. 

We want to apply the gauge-fixing proposal to 
other previous proposals of a LxGT which sup- 
posedly failed due to lack of gauge invariance. For 
this purpose we have chosen the waveguide for- 
mulation of the domain wall fermion where, with- 
out gauge fixing, mirror chiral modes appeared at 
the waveguide boundary in addition to the chiral 
modes at the domain wall or anti-domain wall 
to spoil the chiral nature of the theory |5j. Our 
investigation of the reduced model in this case re- 
veals that both at Yukawa coupling y = 1 ||] and 
y <C 1 the scalars are fully decoupled, there 
are no mirror modes and the spectrum is that of 
free domain wall fermions. 
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2. Gauge- fixed domain wall action 

For Kaplan's free domain wall fermions ^ on a 
4 + l-dimensional lattice of size L*Ls where Ls is 
the 5th dimension, with periodic boundary con- 
ditions in the 5th or s-direction and the domain 
wall mass m(s) taken as 



m{s) — 



-Too 


Too 



0<s<if 
^ <s < Ls 



the model possesses a lefthanded (LH) chiral 
mode bound to the domain wall at s = and a 
righthanded (RH) chiral mode bound to the anti- 
domain wall at s = For moLs 3> 1, these 
modes have exponentially small overlap. 

A 4-dimensional gauge field which is same for 
all s-slices can be coupled to fermions only for 
a restricted number of s-slices around the anti- 
domain wall with a view to coupling only to 
the RH mode at the anti-domain wall. The gauge 
field is thus confined within a waveguide, WG = 
[s : sq < s < si) with sg = "^°4^^ — 1, Si = 
•^-^"^^ — 1. For convenience, the boundaries at 
(so7 sq + 1) and (si, Si -I- 1) are denoted waveguide 
boundary / and II respectively. The symmetries 
of the model remain exactly the same as in Q . 

Obviously, the hopping terms from sq to sq + 1 
and that from si to si 4- 1 would break the lo- 
cal gauge invariance of the action. This is taken 
care of by gauge transforming the action and 
thereby picking up the pure gauge dof or a ra- 
dially frozen scalar field (Stiickelberg field) at 
the WG boundary, leading to the gauge-invariant 
action (with Lp^ —>■ gxfxi where g € gauge group): 



Sf = 



E 

seWG 



{Ip{U)-W{U) + m{s))^'' 



+ E V''* w + to(s)) + ^ V'V'* 



si^WG 



J2 {i^'PLr^' + i^'^'pRr 

- y i^"' (jyPLr'^^ Pr^') (1) 

where we have taken the lattice constant a = 1 
and have suppressed all other indices than s. The 



projector Pl{r) is (1 T 75)/2 and y is the Yukawa 
couphng at the WG boundaries. The p [U) and 
W{U) are the gauge covariant Dirac operator and 
the Wilson term (with Wilson r = 1) respectively. 
^ and w are the free versions of If) {U) and W{U) 
respectively. 

The gauge- fixed pure gauge action for U(l), 
where the ghosts are free and decoupled, is: 

SB{U)^Sg{U)+ Sgf (U) + Set (U) (2) 

where, Sg is the usual Wilson plaquette action; 
the gauge fixing term Sgf and the gauge field 
mass counter term Set are given by (for a dis- 
cussion of relevant counterterms see f|llP|), 



Sct{U) 




u{u)xyn{u\ 




where 0[U) is the covariant lattice laplacian and 
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Bx = 



^nx — n ^" ^nx 



(5) 



with V^x = i {U^x - Ulx) and k = l/{2^g^). 
Sgf has a unique absolute minimum at U^x = 
1, validating weak coupling perturbation theory 
(WCPT) around g = or k = oo and in the naive 
continuum limit it reduces to ^ J d*x(9^^^)^. 

Obviously, the action SsiU) is not gauge in- 
variant. By giving it a gauge transformation 
the resulting action SB{4>\U^x<f)x+ii) is gauge- 
invariant with Ufj^x gxU^xgl+jx and 4> X ^ gxYx^ 
gx G U{1). By restricting to the trivial orbit, we 
arrive at the so-called reduced model action 

Sreduced = SpiU = 1) + Sb{4>1: 1 ^i'+a) (6) 

where Sf{U = 1) is obtained quite easily from 
eq.(P and 

Sb{^11 ipx+i,) = -K.^ipl{Utp)x 



[^iia^p)x-Bl] 



(7) 



now is a higher-derivative scalar field theory ac- 
tion. Bx in (u\) is same as in (|) with 



(8) 
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3. Weak Coupling Perturbation Theory in 
the reduced model 

At y = 1 , we carry out a WCPT in the coupling 
1/k for the fermion propagators to 1-loop. 

3.1. Free propagators 

In order to develop perturbation theory, in re- 
duced model, we expand, (j)x — exp{ib9x:) where 
b = 1/V2k, leading to free propagator for the 
compact scalar 9 M, 



1 



K 

k 



(9) 



where, fc^ = 2sin(fc^/2). 

Free fermion propagators at y = 1 are obtained 
in momentum space for 4-spacetime dimensions 
while staying in the coordinate space for the 5th 
dimension following (results in cannot be 
directly used because of difference in finer details 
of the action). The appropriate free action to 
start from is. 



5" 



p,s.t 



ii>5s^t + MPl + M^Pr 



(10) 



where, M = F{p)5sa+Mq, Mq = a{s)Ss^t~Ss+i^t, 
a{s) = 1 + m(s), F{p) = ^^^(l ^ cos(j3^)), 

= sin(pp) and — ^^p^- The free fermion 
propagator can formally be written as, 

A(p) = + MPl + Mtp^l 

= {-i} + aP)PlG^ + {-ii> + M)PrG^ (11) 
Exphcit solution of G^'-^^p) are obtained from 



[f + l + Bisf] G^^, - Bis . 
^B{s)Gti.t = Ss,t 



1)G 



s+l,t 



(12) 



and a similar equation for G^ . In (12), B{s) — 
F{p) + l+m{s). We show only the calculations for 
obtaining G^ and henceforth drop the superscript 
L. Setting G = in the region < s < Ls/2-1 
where B{s) = F{p) + 1 — rriQ = a_ and G = G+ 
in the region Ls/2 < s < Lg — 1 where B{s) ~ 
F{p) + 1 + mo = a+. 



f 1 + a 



(13) 



{p' + l + <) G+ - a+G+.i,, - a+G+„i,t 

= 5s,t. (14) 

The solutions of these equations are expressed 
as sum of homogeneous and inhomogeneous so- 
lutions: 

cosh[a±(p)(|s-i|-V2)] 



(15) 



2a± sinh(a±(p)) sinh(a±(p)//2) 

where, cosh(a±(p)) = \ {a± + The third 

term in the above is the inhomogeneous solu- 
tion. In order to get the complete solution we 
need to determine the unknown functions g^±{t) 

and g±\t) in ([l5|), which are obtained by con- 
sidering boundary conditions from eqs. ( [l3| ) at 

s = 0, 1, L./2-1, LJ2 + 1, - 1, 



ooGo.tb) =a+Gl^t{p), 

[p^ + l + aD G^Ap)-a^GlM 

= So,t + aoGt _i tip), 



it 



G 



Ls/2-lJ.- 



(16) 



with B{s) = F[p) + 1 = flo at s = 0, 1^/2. Us- 
ing the boundary conditions (^6|) we arrive at an 
equation of the form. 



A.g(t)=X(t), 
where, g(t) = (^i^^ 



(17) 



,(2) 



7^^^) is a 4- 



5V a'l' 

component vector, X(i) is another 4-component 
vector and A is a 4 x 4 matrix. 

The solution to the equations (jl^) is very com- 
plicated in general, particularly for finite L^, how- 
ever g±{t) can be obtained by solving the above 
equations numerically for different t values. This 
way we can easily construct the free chiral prop- 
agators at any given s-slice, including s = t = 
0, Ls/2. The solutions for Gff.{p) and the result- 
ing propagators can be obtained in a similar way. 

3.2. Tree level fermion mass matrix 

Another issue of interest at the free fermion 
level is the spread of the wavefunctions of the two 
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chiral zero mode solutions along the discrete s- 
direction and their possible overlap. A finite over- 
lap would mean an induced Dirac mass. The ex- 
tra dimension can be interpreted as a flavor space 
with one LH chiral fermion, one RH chiral fermion 
and {Ls/2 — 1) heavy fermions on each sector of 
1 < s < Ls/2 - 1 and Ls/2 +1 < s < Ls-1. We 
consider flavor diagonalization of MqMq {Mq = 
M{p = 0) is not hermitian): 



MqM, 



(0)^ 
(0)' 



(18) 
(19) 



where the index j for the eigenvalues and the 
eigenvectors is a flavor index. 

We have carried out explicit solutions for the 
chiral modes at the domain and anti-domain walls 
and for the heavy modes Q. We do not present 
these results explicitly here except to point out 
that the overlap of the opposite chiral modes at 
the domain and anti-domain walls is exponen- 
tially damped anywhere in the s direction. 

3.3. 1-loop fermion self-energy 

Next we calculate the fermion propagators to 
1-loop. A half-circle diagram which is diagonal in 
flavor space contributes to LL and RR propaga- 
tor self-energies and a flavor off-diagonal tadpole 
diagram produces the self-energy for the LR and 
RL parts. However, the self-energies are nonzero 
only at the waveguide boundaries / and //. Be- 
sides, there is also a contribution to the fermion 
self-energy for the LR and RL parts coming from 
a flavor off-diagonal half-circle diagram (we shall 
call it a global-loop diagram) where the scalar field 
goes around the flavor space connecting fermions 
at the waveguide boundaries / and //. 

By expanding (px = exp(i60j:) and retaining 
up to 0{b^) in the interaction term we find the 
vertices necessary to calculate the fermion self- 
energy to 1-loop. 



R 



R ' 



p, So + 1 p- k, So p, So + 1 



The LL propagator on the (sq -I- l)-th slice at 
the waveguide boundary / receives a nonzero self- 
energy contribution from the half-circle diagram, 

-(SL(P)),,= / b'Hj,{q),PLGL{q)l^^^^^ 
J k 

g{k) Ss^s„+iSt,so+i (20) 

^7TEf4«(9)l g{k) 6s,so+i5t,so-,i (21) 

1j L J So, So 

k 

where q = p — k and the expression in the 
square bracket in (|2^) is the free RR propaga- 
tor [5j^^]so,so on the so-slice. Eq.(|20|) assumes 
infinite 4 space-time volume while in eq. ( ^l|) a fi- 
nite space-time volume is considered. To avoid 
the infra-red problem in the scalar propagator, we 
use anti-periodic boundary condition in one of the 
space-time directions in evaluating ([2l|). 

In a similar way, 1-loop corrected RR or LL 
propagators are obtained at all the s-slices of the 
waveguide boundaries / and //, i.e., at the slices 
So, So -I- 1, si and si -I- 1. 

^ ' *" ~ 

/ k \ 



R 



So + 1 



P, So 



For the LR propagator connecting so and so -|- 
1 at the waveguide boundary I, the self-energy 
contribution from the tadpole diagram is given 
by, 



(Si«(p)) 
1 
2 
1 



= ^b^Pi 



1 



k k^{k'^ + w2) 



^s,so^t,so + l (22) 

(23) 



where T ~ 0.04 is the tadpole loop integral. 

Similarly the self-energy contribution to the 
LR propagator at the waveguide boundary // 
connecting si and si -I- 1 comes from a tadpole 
diagram and is given by, 

- (si^K(p))st = \b^PLr 5s^sA.s,+i. (24) 
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The global-loop diagram originates from the 
fact that the tp field that couples the fermions 
at the waveguide boundary / is the same field 
coupling the fermions at the waveguide boundary 



' R 



R 



p, So + 1 'So p - k si + 1 



Self-energy contribution from the global-loop 
diagram for the LR propagator is: 



g{k) 5s,sA.so+i (25) 
b''PLn5s,sA.so + l (26) 



where q — p — k and TZ is the loop integral in 
eq.(§|). _ 
The mass parameter Mq gets modified to Mq 
at 1-loop as: 

iMo)stPL = (Afo).tPL + [-(Sifl(0)).t] 

+ [-{^'Uo))st] + [-(si«(o)).. 

= iMo)stPL+b''i'SLR)stPL- (27) 

^LL,flfl(0) = identically. {M^)stPR gets modi- 
fied accordingly. 

3.4. Fermion mass matrix diagonalization 
in 1-loop 

Diagonalization of the fermion mass matrix 
at 1-loop 1^ shows that, i) the zero modes are 
perturbatively stable, and ii) the overlap of the 
opposite chiral modes at the domain and the 
anti-domain walls are still exponentially damped. 
This clearly rules out the necessity of a fermion 
mass counter term. 

4. Numerical results at y = 1 

In the quenched approximation, we have first 
numerically confirmed the phase diagram in 
of the reduced model in (k, k) plane. The phase 
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Figure 1. Schematic quenched phase diagram. 

diagram shown schematically in Fig.l has the in- 
teresting feature that for large enough k, there 
is a continuous phase transition between the bro- 
ken phases FM and FMD. FMD phase is char- 
acterized by loss of rotational invariance and the 
continuum limit is to be taken from the FM side 
of the transition. In the full theory with gauge 
fields, the gauge symmetry reappears at this tran- 
sition and the gauge boson mass vanishes, but the 
longitudinal gauge do/remain decoupled. In Fig.l 
PM is the symmetric phase and AM is the broken 
anti-ferromagnetic phase. The numerical details 
involved in reconstruction of the phase diagram 
and the fermionic measurements that follow will 
be available in [Q. 

For calculating the fermion propagators, as in 
Q we have chosen the point k = 0.05, k = 0.2 
(gray blob in Fig.l). Although this point is far 
away from k — oo, around which we did our per- 
turbation theory in the previous section, the im- 
portant issue here is to choose a point near the 
FM-FMD transition and away from the FM-PM 
transition. The results below show that for the 
fermion propagators there is excellent agreement 
between numerical results obtained at k = 0.05, 
K = 0.2 and perturbation theory. 

Numerically on 4^^ x 16 and G'^ x 16 lattices with 
Ls = 22 and too = 0.5 we look for chiral modes 
at the domain wall (s = 0), the anti-domain wall 
(s = 11), and at the waveguide boundaries (s = 
5,6 and s = 16, 17). Error bars in all the figures 
are smaller than the symbols. 
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Figure 2. Propagators at domain wall s = 
(Ls = 22; a.p.b.c. in L^, y = 1.0). 



Figure 3. Chiral propagators at domain wall s = 
(L, = 22; a.p.b.c. in L4, y = 1.0). 



Figs. 2 and 3 show the RR propagator |S'_r_r| 
and the LL propagator \Sll\ at the domain and 
anti-domain wall as a function of a component 
of momentum p4 for both p = (0, 0, 0) (physical 
mode) and (0, 0, tt) (first doubler mode) at y = 1. 
From the figures, it is clear that the doubler does 
not exist, only the physical RR (LL) propagator 
seems to have a pole at p = (0, 0, 0, 0) at the anti- 
domain (domain) wall. In all the figures, NS, PT 
and FF respectively indicate data from numerical 
simulation, from perturbation theory and from 
free fermion propagator by direct inversion of the 
free fermion matrix. 

For Figs. 2 and 3, PT also mean zeroth order 
perturbation theory, i.e., numerical solution of 
propagator following eq.(p7|). We have PT re- 
sults also for 6'^ x 16 lattice but have chosen not 
to show them because they fall right on top of 
the numerical data. Instead PT results shown for 
8"^ X 20 lattice for which the p4 points are dis- 
tinct. The dotted line in all figures refer to the 
propagator from PT using a 256'^ x 1024 lattice. 
The curves stay the same irrespective of methods 
or lattice size. Based on the above, we can con- 
clude that there are only free RH fermions at the 
anti-domain wall, and at the domain wall there 
are only free LH fermions. 

Figs. 4 and 5 show no evidence of a chiral mode 
at the waveguide boundaries s = 5 and 6 and ex- 



cellent agreement with 1-loop perturbation the- 
ory. Here too doublers do not exist. Actually the 
agreement with the FF method (direct inversion 
of the free domain wall fermion matrix on a given 
finite lattice) is also excellent, because the 1-loop 
corrections are almost insignificant. For clarity, 
in Figs. 4 and 5 we have not shown the LL prop- 
agator on s = 5 and the RR propagator at s = 6, 
but conclusions are the same. 

Similar investigation at the other waveguide 
boundary s — 16, 17 also does not show any chi- 
ral modes. Previous investigations of the domain 
wall waveguide model without gauge fixing Q 
have shown that the waveguide boundaries are 
the most likely places to have the unwanted mir- 
ror modes. This is why we have mostly concen- 
trated in showing that there are no mirror chi- 
ral modes at these boundaries, although we have 
looked for chiral modes everywhere along the fla- 
vor dimension. In fact, we do not see any evi- 
dence of a chiral mode anywhere other than at 
the domain wall and the anti-domain wall. 

5. Numerical results at j/ <C 1 

At y = 0, the domain wall and the anti-domain 
wall are detached from each other. Fermion 
current considerations and numerical simulations 
clearly show in this case that mirror chiral modes 
form at the waveguide boundaries. 
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Figure 4. i?i?-propagator at waveguide boundary 
s = 5 {Ls = 22, a.p.b.c. in L4, ?/ = 1.0). 

On the other hand, at ?; = 1 (as presented 
above) and y > 1 |l2[ , the mirror chiral modes 
are certainly absent at the waveguide boundaries. 
ActuaUy in this case, the only chiral modes are 
at the domain wall and at the anti-domain wall 
and the spectrum is that of a free domain wall 
fermion. 

The interesting question is what happens at y 
smaller than unity, especially at positive values 
near zero. To investigate this question, we looked 
for chiral modes at y = 0.75, 0.5, 0.25 at places 
other than the domain wall and the anti-domain 
walls, especially at the waveguide boundaries. On 
a 6^ X 16 lattice, the chiral propagators at the 
waveguide boundaries showed an increasing trend 
as the fourth component of momentum p4 (with 
p = 0) was decreased to the minimum value pos- 
sible for this lattice size, something that could 
signify a pole at zero momentum. To resolve this 
we took lattice sizes which were bigger in the 4- 
direction to accommodate lower p4-values. We 
found that for each y < 1 there is a big enough 
lattice size for which the chiral propagators ulti- 
mately start showing a descending trend. Obvi- 
ously the smaller the Yukawa coupling, the big- 
ger the L4 extension was required. Moreover, all 
the chiral propagators at these small Yukawa cou- 
plings matched exactly with the corresponding 
free case. Our results at the waveguide bound- 



Figure 5. iL-propagator at waveguide boundary 
s = 6 {Ls = 22, a.p.b.c. in L4, y = 1.0). 

aries s = 5, 6 are summarized in the Figs. 6 and 7 
for the RR and the LL propagators respectively. 
Poles for these chiral propagators are clearly not 
seen. 

6. Conclusion 

Using the gauge fixing approach and tuning 
only a finite number of counterterms (in this case, 
just the K-term), we end up in the reduced model 
with free domain wall fermion theory consisting 
only of LH and RH chiral modes respectively at 
the domain and the anti-domain wall. With the 
U(l) transverse gauge do/back on the waveguide, 
only the RH fermions on the anti-domain wall will 
be gauged (according to our construction). The 
gauge degrees of freedom are completely decou- 
pled. 

We reached our conclusions for the case of 
y = 1 by performing a perturbation theory for 
the fermion propagators around k = 00 and com- 
paring them to the numerical results at k = 0.2. 
The comparison is near perfect due to the robust 
properties of domain wall fermions (this does not 
happen nearly as nicely for the Wilson fermions) 
and shows that as long as R is big enough to avoid 
the FM-PM transition, it is as good as infinity, 
i.e., the whole FM-FMD transition line is in the 
same universality class as the perturbation point 
K = 00. 
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Figure 

s — 5 with (a) j/ 
0.75, (d) y ^ 1.0 (Ls 22 



(b) y = 0.50, (c) 
a.p.b.c. in L4) 



Investigations at y ^ 1 lead to the same qual- 
itative conclusions and indicate that the model 
for any nonzero Yukawa coupling belongs to one 
universality class. 

The transition to a nonabelian gauge group 
in this gauge-fixing approach is nontrivial and 
should be pursued. A more detailed account of 
our studies can be found in . 
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Figure 7. LL propagator at waveguide boundary 
s = 6 with (a) y = 0.25, (b) y = 0.50, (c) y = 
0.75, (d) y = 1.0 {L, = 22; a.p.b.c. in L^). 
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